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Abstract 

Glauber dynamics is a powerful tool to generate randomized, approximate solutions 
to combinatorially difficult problems. Applications include Markov Chain Monte Carlo 
(MCMC) simulation and distributed scheduling for wireless networks. In this paper, 
we derive bounds on the mixing time of a generalization of Glauber dynamics where 
multiple vertices are allowed to update their states in parallel and the fugacity of each 
vertex can be different. The results can be used to obtain various conditions on the 
system parameters such as fugacities, vertex degrees and update probabilities, under 
which the mixing time grows polynomially in the number of vertices. 

1 Introduction 

Consider a graph G = {V,E), where V is the set of vertices and E is the set of edges. 
Suppose \V\ = n. For each vertex f € F, we use Nv = {w : {v,w) ^ E} to denote the 
set of neighbors of v in the graph. An independent set of G is a subset of the vertices where 
no two vertices are neighbors of each other. Let X be the set of all independent sets of G. 

A configuration of the vertices in G is a vector of the form {ay)y(^v, with cTt, G A = {0, 1} 
for all f € V . For a vertex v and a configuration a G A"', we say w € cr if cj^, = 1. A 
configuration a on G is feasible if the set {w € F : cxt, = 1} is an independent set of G, i.e., 
if 

(^v + (^w < 1, for all {y^w) G E. (1) 

Let Vt C A"- be the set of all feasible configurations on G. 

We associate each vertex v with a parameter A. We are interested in the following 
product-form distribution over the feasible configurations (independent sets) of the graph: 



vdcr 

z = VTTa. (3) 



o-GSI veer 



Note that this corresponds to the so-called hard-core gas model studied in statistical physics, 
where A is called the fugacity (e.g., [6l[9]). 
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Glauber dynamics is a Markov chain which generates the stationary distribution in ([2]). 
It has many applications in statistical physics and computer science (e.g., hard-core gas 
model, graph coloring, approximate counting, combinatorial optimization [2l[3l[6]). Under 
(single-site) Glauber dynamics, in each time slot one vertex is selected uniformly at random, 
and only that vertex can change its state while other vertices keep their states unchanged. 
Let a{t) be the state of the Markov chain in time slot t. 



Single-Site Glauber Dynamics (in Time Slot t) 



1. Choose a vertex v (zV uniformly at random. 

2. For vertex v: 

If E»6A^„^«-(i-l) = 

(a) (Ty{t) = 1 with probability p = 

(b) fT„(t) = with probability p = 
Else 

(c) a^{t) = 0. 

3. For any vertex w &V \ {v}: 

(d) au,{t) =awit-l). 



It is not hard to verify that the Glauber dynamics Markov chain is reversible and 
has the product-form distribution in ([2]). In most applications, the performance of the 
Glauber dynamics is determined by how fast the Markov chain converges to the stationary 
distribution. The Glauber dynamics is said to have the fast (rapid) mixing property if the 
mixing time is polynomial in the size of the graph (the number of vertices n). In [9] it was 
shown that single-site Glauber dynamics has a mixing time of O(nlogn) when A < -^r^, 
where A is the maximum vertex degree in the graph. 

Recently, Glauber dynamics has been applied to design distributed throughput-optimal 
scheduling algorithms for wireless networks (e.g., OEIIS]). In the wireless network setting, 
the graph G = {V, E) corresponds to the interference graph of the wireless network, where 
the vertices in V represent links (transmitter-receiver pairs) in the network, and there is 
an edge between two vertices in G if the corresponding wireless network links interfere 
with each other. A feasible schedule of the network is a set of links which do not interfere 
with each other, which corresponds to an independent set in the interference graph G. To 
achieve maximum throughput, the fugacities need to be chosen as appropriate functions of 
the queue lengths of the links, which are normally different from link to link. This motivates 
the study of Glauber dynamics with heterogenous fugacities. 

In addition, in wireless networks, potentially multiple network links (vertices in G) can 
update their states in a single time slot, and we would expect that the mixing time of the 
Glauber dynamics Markov chain will be reduced with such parallel updates. However, the 
Markov chain may not even have the product-form distribution (which is a key property for 
establishing throughput-optimality in [5l[71[8]) if we let an arbitrary set of vertices update 
their states. The following two questions then arise: 
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(1) How to select the vertices in each time slot to update their states such that the product- 
form distribution is maintained? 

(2) What is the mixing time of such a Glauber dynamics with parallel updates? 

The first question has been addressed in [7] and the second question will be addressed 
in this paper. The rest of the paper is organized as follows. In Section [2] we introduce a 
generalization of Glauber dynamics with parallel updates and heterogenous fugacities. In 
Section [3] we provide some technical background on the mixing time of Markov chains. In 
Sections [4] and El we derive bounds on the mixing time of the Glauber dynamics with parallel 
updates and heterogenous fugacities. The paper is concluded in Section [6l 

2 Glauber dynamics with parallel updates 

In [7] we have introduced a generalization of Glauber dynamics where multiple vertices 
(wireless network links) are allowed to update their states in a single time slot, under which 
the Markov chain is reversible and retains the product-form distribution. The key idea is 
that in every time slot, we select an independent set of vertices m € Z to update their states 
according to a distributed randomized procedure, i.e., we select m G X with probability Qm, 
where X^mGX'?m ~ ^^^^ ^ update set (or decision schedule in [7J). The parallel 

Glauber dynamics is formally described as follows. 



Parallel Glauber Dynamics (in Time Slot t) 



1. Randomly choose an update set m € I with probability qr 

2. For all vertex v € m: 

(a) <Tt,(t) = 1 with probability p^ — 



1+A„ • 

(b) (T„ (t) = with probability = • 
Else 

(c) a^{t) = 0. 

For all vertex w ^ m : 

(d) fT^(t) = aw{t-l). 



The following results on the parallel Glauber dynamics have been established in 

Lemma 1 Let m{t) be the update set selected in time slot t. If CT{t — 1) E and m(t) G T, 
then a{t) G 17. 

Because a{t) only depends on the previous state a{t — 1) and some randomly selected 
update set m(t), a{t) evolves as a discrete-time Markov chain (DTMC). Next we will derive 
the transition probabilities between the states. 

Lemma 2 A state a G 11 can make a transition to a state rj £ 0, if and only if a Li rj Q 
and there exists an update set m G X with gm > such that 

a A rj = {a \r]) U {r]\a) Cm, 
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and in this case the transition probability from a to rj is given by: 

pia,v) = E n p-){ Up-){ n p-){ n p-)- 

(4) 

Under the Glauber dynamics with parahel updates, let 
be the probability of updating vertex f in a time slot. 

Theorem 1 A necessary and sufficient condition for the Markov chain of the parallel 
Glauber dynamics to be irreducible and aperiodic is Umei:q^>o^ = V, or equivalently, 
Qv > for all V € V , and in this case the Markov chain is reversible and has the following 
product-form stationary distribution: 

vEcr 

Z = EII^- (6) 

Remark 1 The single-site Glauber dynamics can be viewed as a special case of the parallel 
Glauber dynamics in which > if and only if the update set m always consists of only 
one vertex. 

In this paper we will show that the parallel Glauber dynamics has a very fast mixing 
time O(logn) under certain conditions for bounded-degree graphs. On the other hand, it 
was shown in [4J that the single-site Glauber dynamics has a mixing time at least r2(n log n) 
for bounded-degree graphs. 

3 Mixing Time of Markov Chains 

Consider a finite-state, irreducible, aperiodic Markov chain (P, Q, vr) where P denotes the 
transition matrix, denotes the state space, and tt denotes the unique stationary distribu- 
tion of the Markov chain. 

Definition 1 The variation distance between two distributions onVt is defined as 

\\ll~v\\var = ]^^\l^{x)-v{x)\. (7) 

Definition 2 The mixing time Tmix{^) for e > of the Markov chain is defined as the time 
required for the Markov chain to get close to the stationary distribution. More precisely, 

Tmix{e) =maxinf it : \\P\x,-) -7r||^ar) < (8) 
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Definition 3 A coupling of the Markov chain is a stochastic process {X{t),Y {t)) on $7 x Q 
such that {X{t)} and {Y{t)} marginally are copies the original Markov chain, and if X{t) = 
Y{t), then X{t + 1) = Y\t + 1). 

Let $ be a distance function {metric) defined on 17 x which satisfies that for any 
x,y, z G $7: 

(1) <I>(x, y) > 0, with equahty if and only if x = y. 

(2) = $(y,x). 

(3) $(x,z) < ^x,y)+^{y,z). 

Let ^ 

Drain = miu ^{x,y), Dmax = max $(x,y), D = -Jn^, 

x,y£n,x^y x,y&il ^Jmin 

The following result can be used to obtain an upper bound for the mixing time of the 
Markov chain (e.g., [2]). 

Theorem 2 Suppose there exist a constant /? < 1 and a coupling {X{t), Y{t)) of the Markov 
chain such that, for all x,y E 

EmX{t + l),Y{t + l))\X{t) = X, Y{t) = y]< /3<^{x, y). (9) 

Then the mixing time of the Markov chain is bounded by: 

TmUe) < i^^^^. (10) 

In general, determining /3 is hard since one needs to check the contraction condition ([9]) 
for all pairs of configurations. In [T] the so-called path coupling method was introduced by 
Bubley and Dyer to simplify the calculation. Under path coupling, we only need to check 
the contraction condition for certain pairs of configurations. The path coupling method is 
described in the following theorem. 

Theorem 3 Let S Q 0, x 0, and suppose for all X,Y x 0,, there exists a path X = 
Zq, Zi, . . . , Zr = Y between X and Y such that {Zi, Zi^i) € S for < I < r and 



<^{X,Y) = Y,^Zi,Zi+i] 



1=0 

Suppose there exist a constant /? < 1 and a coupling {X(t),Y{t)) of the Markov chain such 
that for any (x, y) G S, 

E[^{X{t + 1), Y{t + l))\X{t) = X, Y{t) =y]< /3c^(x, y). 

Then the mixing time of the Markov chain is bounded by: 

r„,.(.) < (11) 

Note that the key simplification in the path coupling theorem (Theorem [3D , compared 
to the coupling theorem (Theorem [2]) , is that the contraction condition ([9]) needs to hold 
only for (x, y) € S, instead of (x, y) G Q x Vt. 
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4 Mixing Time of Glauber Dynamics With Parallel Updates 



In this section we analyze the mixing time of Glauber dynamics with parallel updates using 
the path coupling theorem. We will use the following distance function: for any o", G 0, 

$((7,7?) = Y,Wv-Vv\f{v)= /W' (12) 

where f{v) > is a (weight) function of v (z V and recall that a A rj = (a \r]) L) {r]\a). 
Note that this distance function is a weighted Hamming distance function and satisfies all 
the properties of a metric. 

Consider the following coupling {a{t),r]{t)): in every time slot both chains select the 
same update set and use the same coin toss for every vertex in the update set if that vertex 
can be added to both configurations. 

Let E[A^{a{t),r]{t))] be the (conditional) expected change of the the distance between 
the states of the two Markov chains {cr(i)} and {r]{t)} after one slot: 

E[A<i>iait),r^m = EMait + l),r/(t + l)|(7(t), r?(t)] - $(a(t), ??(t)). 

For any m € X, let 

E[A"'^{a{t),r]{t))] = E[A^{a{t),r]{t))\m is the update set]. 

Lemma 3 Let lii = (yi, . . . ,y\m\) '^'^y ordering of m. For any a{t),r]{t) € 0, 

m 

E[A^^{a{t),^m = Y,m"'^i^it)Mm. 

k=l 

Proof Note that the value of $((T, r/) is completely determined by the set a Arj, which in 
turn depends only on a and r]. Hence, it suffices to show that we will obtain the same sets 
a{t + 1) and rj{t + 1) by updating all the vertices of m simultaneously and by updating 
them in any sequential order. 

The moves trying to remove vertices from a configuration will be successful in all cases. 
The outcome of a move trying to add a vertex y to a configuration u is successful if and 
only if Afy H a; = 0. But m is an independent set, so no neighbor of a vertex y G m is 
in m. Then the states of the neighbors of y are unchanged after updating any subset of 
vertices of m. Hence, V/i;G{l,...,|m|}, the outcome of a move trying to add y^ to uj will 
be the same if we update the vertices of m sequentially and if we update all the vertices 
of m simultaneously. This is in particular true for the configurations a{t) and T]{t), so we 
conclude that the two update procedures will yield the same sets a{t + 1) and rj{t + 1). 

We say that a,r] £ Q, are adjacent and we write a ~ 77 if there exists v G V such that a 
and rj differ only at v. Let 

= I (fj, 77) : (T, 7? € $7 and a r]^. 

Note that under the distance function defined in (I12p . for all a,r] € 17, we can find a 
path fj = To, Ti, . . . , T\crAr]\ = V between a and rj such that (t/, r^+i) S S for < / < r and 

Now consider a pair of adjacent configurations a{t) and 7?(t) that differ only at v. With- 
out loss of generality, suppose cj„(t) = and r]v{t) = 1. This means that, r]w{t) = for all 
w G Afy Since, a{t) and r]{t) differ only at v, this also means that aw{t) = for all w G Afv 
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Lemma 4 



E[A^ait),vm < -qj{v) + TT^fi^)- (13) 



Proof Using Lemma [3l we have 

= Y.^rnE[A^^a{t),vm 

m 

m J/Gm 

= 5;(?,ii;[A^<J>(a(t),r?(t))]. 

Note that only updates on vertices v and ii; G A/",; can affect the value of E[A^{a{t),rj(t))]. 
In particular, if v is selected for update and since we use the same coin toss for both Markov 
chains, cr{t + 1) = r/(t + 1). Thus E[A''<^{a{t),rj{t))] = -f{v). 

If w G Mv is selected for update, under configuration r/(t), w can only take value 
because w has a neighbor (i.e., v) belongs to r]{t). While under configuration a{t), there 
are two cases: 

1) if w has a neighbor in cr(t), then w can only take value 0; 

2) if w has no neighbors in cr{t), w can take value 1 with probability j^j;- and value 
otherwise. 

Hence for w G Mv, 

i?[A-cI>(a(t),r?(t))]<-^/(^). 
Summing up all contributions we have (|13p . 



Now we are ready to present the main result of this paper. 
Theorem 4 For any positive function f{v) ofv G V, letm = min^gy f{v), M = maxy^v f{v), 



andi = ^. If 



e^minL^fiv)- ^f{w)\ > 0, (14) 

then the mixing time of the parallel Glauber dynamics is bounded by: 

M 

Trm.ie) < —log{e-'nO. (15) 

Proof For any pair of adjacent configurations {a{t),ri(t)) G S that differ at some vertex 
V & V, from (fT3l) and (fT4|) we have: 



E[A^a{t),vm <-0< --<^{a{t),v{t)), 
where we use the fact that 

<^{a{t),rj{t))=f{v)<M. 
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Therefore, 



Ema{t + l),rj{t + ma{t),r,{t)] < {l - -)^a{t),7j{t)). 



Then, by applying the path couphng theorem where we let /? = 1 — and D = n^, we 
prove the bound in (fT5l) . 



4.1 Conditions for Fast Mixing of Parallel Glauber Dynamics 

We can now specify the (weight) function / to obtain different conditions on the fugacities 
At,'s for fast mixing. We will show three such examples. 

Corollary 1 Let m = min^.gv' M = max^,gy and ^ = ^- // 

= mill i 1 + A„ - ^ A.„ I > 0, (16) 



then we have 



Trmxie) < ^log(e-ine). (17) 



Proof Choose f{v) = yv€V 



Corollary 2 Let qmin = min^jev- q^, Qmax = max„gy Qy, and = If 



_K 
1 



6 ^ max y — ^ < 1, (18) 



then we have 



, ^ log (e , ^ 

T™.(e) < ■ (19) 



Proof Choose /(v) = ^, V-y € 



Remark 2 Ifq^yc for some constant c > which is independent of the size of the network 
n, then the mixing time is O(logn). In particular, for a hounded- degree graph G where 6 
and A are the minimum and maximum vertex degrees, if = where dy is the degree 
of V, then we have < qv ^ o-iT'd 

-tmix e < ;- log — e n . 

Corollary 3 If < -j^^i fof o,ll v gV, then Tmixi^-) ^ ^^og{e~^n^), for some constants 
M, C and9>0. 

Proof The proof is very similar to that of Theorem HI We can choose f{v) = ^ and let 



M = maxt,ey ^ and 



^ _ max^gy ^ 
min^gy ^ ■ 
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The parallel Glauber dynamics will have fast mixing if 



mm < 



To achieve that we need, for all v (zV, 



dy - "\ d,„ > > 0. 



X,, 



d,,, > 0. 



It is sufficient that j^j^d^ < 1, which is equivalent to A^, < 

Note that the condition A^ < j^rj foi" all v S ^ might be very different from b = 
maxj^gy Z]u;eA/'„ T+X^ ^ (^'S-' ^ star network). 

5 Mixing Time of Single-Site Glauber Dynamics With Het- 
erogenous Fugacities 

Since the single-site Glauber dynamics is a special case of the parallel Glauber dynamics 
(in which the update set always consists of only one vertex), the general results derived in 
the previous section also apply to the single-site Glauber dynamics. The motivation of this 
section is to derive a larger region on the fugacities under which the single-site Glauber 
dynamics is fast mixing, where we use a similar path coupling technique as in [9]. 

We redefine the state space to be the set of all configurations A". For a vertex v 
and a configuration o" € A", we define the set of blocked neighbors of v with respect to a as 

Bcj{v) = G A/'„ : -u; G (7 or TVu, n <T / 0| 

and of unblocked neighbors of v as 

B,{v)=K\B„{v). 

Note that the unblocked neighbors of v are the neighbors of v that can be added to the 
configuration in the next move. If ?; ^ o", we write a'" to denote the configuration that 
differs from a only at v. 

We say that (T,r] ^ A" are adjacent and we write o" ~ r/ if there exists v such that 
a and rj differ only at v. Let 

S = {(a,r/) : 0- ~ 7?}. 

Then, for all a,r] & A", we can define a path in A" between a and ry as a sequence 
{TQ,...,Tr) C A" such that for < i < r,rj ~ Xj+i and tq = o", = r/. Let V{a,r]) be 
the set of all paths in A" between a and ij. 

We will use the following distance function on A" x A": Vo", 77 G A", let 



r-l 



$(a,r/) = min(^^^..._^^)gp(^_^) ^/(ri,ri+i). 



i=0 
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where 



/(c7,cT-) = l + i Yl (20) 

can be viewed as the length of edge {a, a'") according to the adjacent relationship defined 
by S. $ is clearly symmetric, non-negative, zero only when the configurations are identical, 
and satisfies the triangle inequality (because V{a,r]) C V{cr,fi) x V{fj.,r])). Hence, it is 
indeed a metric on A" x A". 

For all v,u G V, we let T{v, u) = Mv riMu be the set of vertices that form triangles with 
V and u. To each vertex v gV, we associate a fugacity A^,. 

Consider a pair of adjacent configurations a (where v ^ a). Suppose a{t) = a, 

r]{t) = a'" . The coupling is simply that each configuration attempts the same move at every 
time slot. More precisely, both Markov chains select the same vertex to update, and use 
the same coin toss if the vertex can be added to the configurations. For convenience, we 
use the following notations from [9]. Let 

E[A<i>] = E[^a{t + l),r,{t + l)\a{t),r,{t)] - cl>(a(t), r?(t)). 

This can be further calculated via the analysis of individual moves. Let 

£'[A+'^<I>] = ii;[A<i>|both chains attempt to add y at time t], 
£^[A^^<1>] = £'[A$|both chains attempt to remove y at time t], 

and denote the total effect of all moves on y by 

E[Ay^] = -^^E[A+y^] + —^E[A-y<^]. 

1 + Ay 1 + Aj^ 

We provide the main result of this section in the following theorem, which says that the 
mixing time of (single-site) Glauber dynamics with heterogenous fugacities is 0(n log n) 
under certain conditions. 

Theorem 5 Let 

a = max < > A^ > (21) 
vGV L ^ — J 

and 7 = "^y^yi^ + A^). If a < 2 and the probability of selecting vertex y & V to update is 
1y ~ then the mixing time of the Glauber dynamics is bounded by 



Trm.{e) < Y^^ogl^e-'n{l + -)j. (22) 

We need the following lemma to prove Theorem^ and its proof is given in the Appendix. 
Lemma 5 

2 ^(1 + Xy)E[Ay<^] < -2 + ^ A^. (23) 
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Proof (Theorem E]) Suppose a < 2 and qy = then using Lemma Owe have 

E[A<^>] = "^qyElAy^] 

yev 

E,6y(l + A,)ii;[A^$] 



< 



-1 + - 



7 
7 

where we use the fact that <I>(cr(t), r/(t)) > 1 for any a{t) ^ rj{t). Therefore, 

EMait + l),r/(f + < (1 + Zi±l)$(a(f),,?(t)). 

7 

Since $(cr, cr'') > 1, so for all a, G A", a ^ i], we have <^(cr, r?) > 1. Also, $(o-, cr") < 
1 + |, and thus ^{a,r]) < n{l + |). Then, apply the path coupling theorem where we let 

13 = 1 + and D = n{l + f ): 

Tm^xie) < ^ log (n(l + g)e-^ 



< ^log(2e-in). 



Example 1 Consider a graph G with n vertices, labelled 1,2, ... ,n. For each vertex i, let 
■^i = {i • J 7^ ^ o,nd \ j — i\ < 3} he the set of neighbors of i. For Glauber dynamics with 
homogeneous fugacity X, Vigoda's result ^9] says that it has mixing time O(nlogn) when 
A < -^7:2 =0.5, since the maximum vertex degree A = 6 in this graph. While for Glauber 
dynamics with heterogenous fugacities, our result says that even the fugacity at some vertex 
exceeds 0.5 (hut is less than 2), the Glauber dynamics can still have 0(n log n) mixing time, 
e.g., when Aj = 1 for vertices i = 6/ + 1, / = 0, 1, . . . , and Xj = 0.18 for other vertices j. 



6 Conclusion 

In this paper we have analyzed the mixing time of a generalization of Glauber dynamics 
with parallel updates and heterogeneous fugacities. By applying the path coupling theorem 
and by choosing appropriate distance functions, we have obtained various conditions on the 
system parameters such as fugacities, vertex degrees and update probabilities under which 
the parallel Glauber dynamics is fast mixing. In particular, we have shown that the mixing 
time of the parallel Glauber dynamics grows as O(logn) in the number of vertices n for 
bounded-degree graphs when the fugacities satisfy certain conditions. 



A Proof of Lemma [5] 

Proof Note that when a < 2, the edge length defined in (pO|) satisfies 1 < l{a,a^) < 2. 
Hence the length of a path of two edges or more is at least 2. This implies that the distance 
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between two adjacent configurations is simply the length of the direct edge between them, 
i.e., 

<^{a,a-) = l{a,cj-) = l + \ ^ A^. 

Only updates on vertices within a distance of 2 hops from v can influence the value of 
$((7, cr"). Hence, 

2^(1 + A,)E[A^$] 
yev 

= 2(1 + A^)S[A''$] + 2 ^ (1 + A^)S[A^$] + 2 ^ (1 + A^)£;[A"$]. 



w 



We will examine separately the terms involving w's and x's. It is important to no- 
tice that the moves trying to remove a vertex from a configuration are always successful, 
whereas the moves trying to add a vertex y to a configuration a succeed only if My H a = 0. 

E[A^*]: 

Only if Mv n cr = 0, will the move trying to add v to cr be successful. Therefore, 



2(1 + A^)E[A''$] = 



-{l + Xv){2 + T.n,&B.{v)^w) ifAA„nc7 = 

-(2 + HweB^iv) ^w) otherwise. 



E[A^*]: 

Trying to update a neighbor of v can only increase the distance between the two config- 
urations. Indeed, either adding a neighbor of v will fail, or it will increase the length of 
any path between the two configurations by one, because it is then necessary to remove 
that neighbor from both configurations before it may be possible to add v. Removing a 
neighbor of v from the configuration will potentially unblock some other neighbors of v and 
thus increase the distance between the configurations. 

More precisely: if if € Bfj{v), it is impossible to add w to a and so £?[A = 0; 
otherwise, £'[A+"'$] = ^{a^j^a'") — ^{a,a^). By the triangle inequality. 

We have ^{aw,CTyj,v) < ^(o", o"'') because Bcr^{v) C Ba{v). Also, 2^{ayj,v,(^'') = 2 + 
^yGB^v (w) K ^"^^ -^CT" (^) = ^'^ ('^) \ {^} \ '^('"' • Then, by combining everything, we get 



2£;[A+"'$] < 



iiweB„{v), 
2 + Eyes, {w)\{v}\T{v,w) \ otherwise. 



Trying to remove w will succeed and unblock some of v's neighbors, as well as add w 
to B„{v) if Nw n (T = 0. A neighbor w' of v will be unblocked if and only if w is the only 
neighbor of w' in a. Thus, 



if t(; ^ fj, 

5{M^ n a = 0)A^ + Y.w'eT{v,w):M^,na={w} otherwise. 
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where we use the indicator function 5 {a) which is equal to 1 if a is true and otherwise. 



As in the case of the terms involving w^s, trying to remove a vertex x from the configurations 
unblocks some neighbors of v and thus increases the distance between the configurations. 
On the other hand, adding a vertex x blocks some neighbors of v and reduces the distance 
between the configurations. More precisely, when removing x from a and a'" , a neighbor w 
of V will be unblocked if and only if x is the only neighbor of w in a. Thus, 

I EweT{v,xy.M^na={x} otherwise. 

Trying to add x fails if x € cr or J\fx fl a 7^ 0. Otherwise, it succeeds and blocks some 
neighbors of v. More precisely, a neighbor ti; of w is blocked if it is not in a, it was not 
already blocked in configuration a before the move, and it is a neighbor of x. Thus, 



2^fA+^$l = I ° if X G a or 

I -T.n^eBM:xeM^^n, Otherwise. 



if X G cj or TVa; n 0- / 0, 

veB^{v):x£Afw 

We can rewrite this equation in the following form: 



w£Bcy{v):x£Ba{w) 



We are now ready to sum up the contributions from all the updates: 
< -(l + A„(5(AA^n(j = 0)) I 2+ Y 

\ w&B„{v) 

+ 5] A. (2+ Yl ^y] (24) 

w&Ba{v) \ y<^Ba(w)\{v}\T{v,w) I 

+ Y ^^^^ ^ = ^^^^ + (25) 

'"Go- Y w' (^T{y,w):M^tr\(y={w} j 

+ E( E (26) 

- E ^4 E ^A- (27) 

x:d{y,x)=2 \weBrj(v):xeBa{w) j 



We want to simplify the terms (i24|) + (j27l) : 



24D + (I2T]) = ^A^ 2+ Y \- ^ 

weB„(v) \ yeB„(w)\{v}\r{v,w) xeB„{w) 
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As we recall that the notation x in {x € Bfj(w)} refers to vertices that are at distance 2 
from V (as opposed to y that refers to any vertex in V), we get that 

{x G B^iw)} = B^{w) \ {v} \ T{v, w), 

that is the set Ba[w) minus its elements at a distance or 1 from v. Then, we can simplify: 

(IMl) + (I27D = 2 ^ A^. 



We now turn to (l25]l + (|26D : 
(ESI) + (IMl) 



= d{Mu, na = 0)A^ + + E E 

= Y^ + Att. < At„. 

Combining all these inequalities, we have 

yev 

< -{l + X^5{Kna = (/})) (2+ A^) + 2 A» + 

= -2 + ^ A^ - XJ{Mv n (T = 0) (^2 + Y 

weAfv weBa{v) 

-2+ Y 



< 
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